G3 Genes|Genomes|Genetics

Laboratory cultivation subjects model organisms to selective pressure and genetic drift that can result in the accumulation of many genomic and phenotypic differences over time. The nematode Caenorhabditis elegans has been used for research since the 1970s, and studies comparing the N2 Bristol and CB4856 Hawaiian isolates provided foundational knowledge about metazoan genome evolution. Most comparative genomics studies have used these isolates because their long-term geographical isolation promoted a high degree of genomic divergence within the species. Further, there is growing evidence of phenotypic differences between laboratory lineages of each wild type isolate after repeated independent lab cultivation of these strains. To examine the genomic divergence between different laboratory lineages the Bristol and Hawaiian backgrounds, we first generated de novo genome assemblies of two Bristol and two Hawaiian lineages from Illumina and PacBio sequencing reads. Following genome assembly, we quantified Single Nucleotide Polymorphisms (SNPs), short insertion/deletions (indels), and genomic structural variants (SVs). Between laboratory lineages of the Bristol isolate, we identified 25,432 SNPs, 5,202 indels, and 441 SVs. When aligning laboratory lineages of the Hawaiian isolate, we identified 4,518 SNPs, 1,188 indels, and 387 SVs. For both sets of comparisons, we find that SNPs and indels are broadly enriched in introns and depleted from coding sequences. In contrast to SNPs and indels, we find that genomic SVs are enriched in intergenic sequences. Taken together, our analyses reveal the accumulation of genomic divergence between lineages of Bristol and Hawaiian C. elegans from independent lab cultivation, and how these variants may underpin emergent phenotypic differences observed in the two most popularly used C. elegans wild type isolates. Author SummaryLaboratory model organisms, like natural populations, are subject to evolutionary pressures and genomic changes during prolonged laboratory cultivation. In this study we comprehensively quantify SNPs, indels, and SVs between independent lab cultivations of the C. elegans lineages of the Bristol and Hawaiian isolates.

The field of genetics of bird migration advances, driven by exponential refinements of sequencing and tracking technologies. In willow warblers (Phylloscopus trochilus), a complex repeat-rich region named MARB (Migration Associated Repeat Block) has recently been found to correlate with the routes taken by individual birds from Europe to their African wintering grounds. However, the genomic location of this region remains unknown. Here, we characterized MARB using a combination of approaches to understand how it evolved. We describe the region using long-read genome assemblies of two willow warbler subspecies (P. t. trochilus and P. t. acredula), two related species, the common chiffchaff (P. collybita) and the greenish warbler (P. trochiloides), and whole genome sequencing data from 76 willow warblers. Finally, we applied karyotyping and fluorescent in situ hybridization techniques on willow warbler spermatocytes to cytogenetically locate MARB. Due to the many repeats, we cannot order scaffolds in silico, but probe hybridization on the karyotype shows that MARB constitutes a single locus (~27.5 Mb) spanning most of the 11th largest chromosome in the willow warbler genome. Interestingly, the MARB regions of all species share several characteristics such as relatively high GC content (50%), a high density of specific repeat families and notably, more than 800 olfactory receptor sequences. Regions homologous to MARB may exist in several migrant bird genomes, though currently unassembled due to their complexity. Resolving these in species with similar migratory polymorphisms to willow warblers will be essential to determine whether MARB influences migratory behaviour across species.

Many taxa have evolved heteromorphic sex chromosomes like the XY system found in mammals. In additional to the sex determination gene which directs development of the gonad into an ovary or testis, sex chromosomes can have drastically different gene content, leading to substantial genetic differences between genetic males and females beyond their gonad identity. Studying the effects of these genetic differences is challenging, as the sex chromosomes and sex determination gene are inherited together, so the effects of genetic differences between the X and Y cannot be easily isolated from the hormonal differences produced by the ovary and testis. The threespine stickleback fish has a heteromorphic XY sex chromosome system and a wide range of well documented sex differences in morphology and behaviors, including complex mating behaviors and male-only parental care. Through genetic manipulation of amhy, the newly identified sex determination gene in threespine stickleback, we are able to generate gonadal males and females with either the XX or XY sex chromosome complement and analyze the separate effects of gonadal sex and sex chromosome complement on sexually dimorphic gene expression. We find that sex chromosomes have a larger effect on gene expression than gonadal sex in somatic tissues, while gonadal sex has a larger effect on expression in the gonads. We also find that the X and Y chromosomes are enriched for genes which show differential expression between females and males. Our findings demonstrate the significant biological impact of sex chromosomes outside of primary sex determination and showcase the utility of the threespine stickleback for studying the genetic basis of sex differences.

At a mere 20 cells, the Caenorhabditis elegans intestine regulates metabolism, energy homeostasis, host defense, yolk production, and genetic aging, all while dynamically responding to its environment. How the intestine develops to carry out these disparate functions is unknown, and how cells differ along the length of the intestine is unclear. To address these questions, we performed single-cell RNA sequencing (scRNA-seq) on FACS-enriched intestinal cells from mixed-stage C. elegans embryos. The resulting single-cell transcriptomes of 974 cells organized into 13 clusters, suggesting a diversity of cell types and states. We used two post hoc approaches to ascribe identities to each cluster. First, genes with known developmental timing in early-, mid-, and late-stages were used to place clusters in time, and smiFISH microscopy was used to fine-tune the assignments. Second, the eight late-stage clusters were assessed for their region of origin. To assign these clusters to anatomical regions, we identified marker genes for each cluster and assessed their expression along the anterior-to-posterior length of the intestine using smiFISH microscopy. Genes associated with growth and cell division were expressed in early stages, whereas genes associated with immune responses and metabolism were expressed later. Genes associated with biotic responses and RNA metabolism were the most likely to vary across the intestines anterior-posterior axis. Finally, perturbation of anterior-localized intestinal transcripts more robustly affected intestinal function compared to central or posterior-localized genes. Overall, this research illustrates the intrinsic heterogeneity across the 20 cells of the embryonic intestine and sets the stage for future works aimed at understanding cell-specific intestinal responses to diet and the environment. ARTICLE SUMMARYWe investigate how the Caenorhabditis elegans intestine develops specialized functions on a spatiotemporal scale. We used single-cell RNA-sequencing to analyze embryonic intestinal cells and identify 13 distinct clusters. Combining gene expression analysis with microscopy, we assigned clusters to developmental stages and anatomical regions. Clusters associated with early intestine development express genes linked to growth and cell division, while later-stage clusters express genes involved in metabolism and immune responses. Genes varied across the intestines anterior-to-posterior axis, and disrupting anterior-specific genes produced stronger functional effects. These findings reveal previously unrecognized intestinal diversity and provide insight into how intestinal cells specialize during development.

Ribosome biogenesis is a highly coordinated pathway that involves the assembly of ribosomal RNAs (rRNAs) with ribosomal proteins (r-proteins) to generate functional ribosomal subunits (r-subunits). The Saccharomyces cerevisiae (yeast) large 60S r-subunit consists of three rRNA molecules and 46 r-proteins. The contributions of nearly all r-proteins of the yeast large r-subunit have been characterized; however, a few non-essential proteins remain poorly understood. Although non-essential, human eL22 has been identified as a key player in p53 regulation during ribosomal stress and as a highly mutated target in cancers. Despite this function, the role of eL22 in ribosome maturation is still ill-defined. In this study, we characterized yeast eL22 r-protein. Our results show that eL22 assembles into intermediate nucleolar pre-60S ribosomal particles. Loss of eL22 impairs cell growth and reduces 60S r-subunit accumulation, phenotypes that are exacerbated at low temperatures. Analysis of pre-rRNA processing by pulse-chase labeling, northern blot hybridization, and primer extension reveals a defect in 27S pre-rRNA maturation, specifically at the level of 27SB pre-rRNA processing. Consequently, nuclear export of eL22-deficient pre-60S particles is mildly impaired. Furthermore, we identify genetic interactions between eL22 and neighboring r-proteins, eL38 and eL31. We conclude that eL22 assembly is required for optimal pre-60S maturation during middle nucleolar stages, particularly at low temperatures, a function likely supported by the cooperative action of other r-proteins associated with common elements of 25S rRNA. HighlightsO_LIWe have studied the role of r-protein eL22 in yeast ribosome assembly. C_LIO_LIeL22 is required for 60S ribosomal subunit production. C_LIO_LIThe absence of eL22 is critical at low temperatures. C_LIO_LIeL22 is important for 27SB pre-rRNA processing and nuclear export of pre-ribosomes. C_LIO_LIeL22 functionally interacts with r-proteins eL38 and eL31 in domain III of 25S rRNA. C_LI

Transcription factor IIIC (TFIIIC) is a multi-subunit protein complex that recruits RNA polymerase III (Pol III) to the majority of its target genes. Evolutionarily conserved overlap between TFIIIC binding sites and the structural maintenance of chromosomes (SMC) complexes suggested a role for TFIIIC in SMC regulation and 3D organization of eukaryotic genomes; but the evidence has remained largely correlational due to the essential role of TFIIIC in RNA Pol III transcription. Here, we directly tested the function of TFIIIC in SMC complex regulation by using auxin inducible depletion in C. elegans. We performed Hi-C and ChIP-seq analyses upon acute depletion of TFTC-3, an essential TFIIIC subunit, and RPC-1, the catalytic subunit of RNA Pol III. Our results show that TFIIIC regulates the localization of two different types of SMC complexes, cohesin and condensin. TFIIIC is also required for increased 3D contacts between distant tRNA genes located on the same or different chromosomes. Depletion of individual SMC complexes did not significantly reduce the intrachromosomal tRNA gene contacts, suggesting redundancy or an independent mechanism mediating these contacts. Together, our study demonstrates an RNA Pol III independent function for TFIIIC, regulating binding of both cohesin and condensin, as well as the 3D organization of tRNA genes.

Current optimum contribution selection (OCS) implementations use point estimates of estimated breeding values (EBVs), potentially leading to suboptimal selections when individuals have uncertain genetic evaluations. We developed a framework assessing how EBV uncertainty affects OCS decisions through MCMC-based approaches using the COSMO optimizer in Julia, evaluated on Norway spruce (Picea abies, n=5,525) and Loblolly pine (Pinus taeda, n=926) populations. Agreement between point estimate (MAP-OCS) and MCMC-OCS was surprisingly low: mean overlap of only 26.6 (4.8) individuals in Norway spruce genotyped subpopulation and 14.1 (3.6) in full pedigree, with Loblolly pine intermediate at 16.0 (9.6). Despite this low individual-level agreement, selection frequency across MCMC iterations corresponded well with EBV rankings (Spearman{rho} = 0.782 for Norway spruce), confirming that higher-EBV individuals were preferentially selected under posterior uncertainty. To comprehensively quantify uncertainty impacts, we employed two complementary metrics: individual robustness scores measuring genetic gain stability upon candidate removal, and population-level contribution distribution metrics capturing concentration of genetic gain across selected individuals. Applying these metrics identified 25 high-risk individuals in Norway spruce and nine in Loblolly pine, and constrained exclusion of these individuals improved individual robustness by 16.5% in Loblolly pine (3.00% genetic gain loss) and 29.8% in Norway spruce (2.14% genetic gain loss). Our uncertainty-aware OCS framework successfully identifies unstable selections that may compromise long-term genetic gain, and we recommend assessing EBV uncertainty through posterior distributions and evaluating population-specific trade-offs when implementing uncertainty-aware selection strategies.

Dauer larvae are a dormant developmental stage in nematodes that is induced by a range of environmental cues. The molecular mechanisms that transduce these cues to regulate dauer entry have been well characterized in Caenorhabditis elegans, whereas those in other nematode species remain unclear. The closest known sibling species of C. elegans, Caenorhabditis inopinata, occupies a distinct ecological niche and shows an extremely low frequency of dauer formation by starvation in laboratory conditions, suggesting that it could serve as a useful comparative model for analyzing dauer-inducing mechanisms. To support such analysis, we generated a fluorescent dauer reporter, Cin-col-183p::mCherry, in C. inopinata based on a previously reported dauer-specific reporter in C. elegans. This reporter showed fluorescence specifically in the pre-dauer and dauer stages, but not in other developmental stages, indicating that it functions as a dauer-specific marker in C. inopinata. Using these marker strains, we compared the responses to high temperature and RNAi-mediated knockdown of insulin/IGF-1 pathway genes (daf-2, age-1, and pdk-1), and found that dauer induction differs mechanistically between C. elegans and C. inopinata. This dauer-specific fluorescent strain will be a useful tool for investigating the diversity of dauer-inducing mechanisms across nematode species. Article SummaryDauer is a dormant developmental stage in nematodes induced by environmental stress. Although its regulation is well studied in Caenorhabditis elegans, the mechanisms in other species remain unclear. Here, we developed a fluorescent dauer reporter, Cin-col-183p::mCherry, in Caenorhabditis inopinata, a close relative of C. elegans. The reporter was specifically expressed in pre-dauer and dauer stages, confirming its usefulness as a dauer marker. Using this strain, we found that responses to high temperature and insulin/IGF-1 pathway gene knockdown differ between C. elegans and C. inopinata. This reporter will help reveal diversity in dauer-inducing mechanisms across nematode species.

Gene model for the ortholog of Lst8 (Lst8) in the May 2011 (WUGSC dyak_caf1/DyakCAF1) Genome Assembly (GenBank Accession: GCA_000005975.1) of Drosophila yakuba. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Interest in quantifying linkage disequilibrium (LD, non-random associations of alleles at different loci) has skyrocketed in recent years as researchers have focused on use of LD in genome-wide association studies (GWAS), for studying historical demography, and for estimating effective population size (Ne). The most widely used LD metric is r2 = the squared correlation of alleles at a pair of loci. Despite a half century of efforts, developing an unbiased expectation of r2 as a function of the many factors that can affect it (physical linkage, genetic drift, selection, migration, mutation, mating systems) remains elusive. Furthermore, even when all of these other factors are absent, empirical estimates of r2 are upwardly biased by sampling a finite number (S) of individuals, and that must be accounted for if one wants to focus on the desired signal of LD. Previous approaches to estimate [Formula] have been shown to be biased to greater or lesser degrees. The purpose of this short paper is to demonstrate that a simple and apparently exact expression for [Formula] does exist for the special case where sampling error is the only factor contributing to r2, in which case [Formula] = 1/(S - 1). When other factors contribute heavily to LD, [Formula] shrinks toward 0 as empirical r2 [->] 1. However, for estimating contemporary Ne with unlinked markers, empirical r2 will generally be small and 1/(S - 1) will provide a robust estimate of [Formula].

Several recent studies have quantified signed linkage disequilibrium (LD) among mutations in genomic datasets, often reporting positive LD, particularly among mutations presumed to be less deleterious, such as synonymous variants. In this article, we investigate two potential sources of this positive LD: the focus on rare alleles, as adopted in several previous studies, and errors arising in the mapping of short-read sequences onto a reference genome. Using coalescent simulations, we extend previous theoretical results of the effect of focusing on rare alleles, and show that derived alleles present at similar frequencies tend to be in positive LD. Reanalyzing datasets from Capsella grandiflora and Drosophila melanogaster, we show that LD among synonymous derived alleles vanishes in the absence of any conditioning on frequency, while LD between mutations categorized as potentially deleterious by the SIFT4G program stays positive. However, we show that in both cases, this positive LD may be at least partly caused by the potential mismapping of a small fraction of sequences in some individuals, which could be a consequence of structural variants that are absent from the reference genome. Overall, these results show that average signed LD among mutations can be strongly affected by technical artifacts even if these concern only a minority of variants. Finally, we discuss other possible sources of positive LD among deleterious mutations.

Animals sleep more when they are sick. In C. elegans, stress-induced sleep (SIS) follows cellular injury such as exposure to ultraviolet (UV) light. The genetic regulators of SIS remain incompletely defined. Using a worm-picking robot, multi-well WorMotel imaging, and association analysis we performed a semi-automated screen of 941 whole-genome-sequenced Million Mutation Project (MMP) strains. We quantified behavioral activity and quiescence before and after ultraviolet (UV) radiation. We applied the Sequence Kernel Association Test (SKAT) to this behavioral data to prioritize 6,663 genes and observed significant enrichment of known SIS genetic regulators. Based on these results, we conducted a candidate validation screen for additional genes regulating SIS. We identified three genes (strd-1, egl-8, cla-1), mutations in which reproducibly influence SIS. Further exploration of these genes holds potential for enhancing our understanding of the molecular basis of SIS. These findings establish a pipeline for automated behavioral phenotyping coupled with gene-based association to accelerate studies of C. elegans neurogenetics.

An expectation of a hypothesis that proposes cell-to-cell signalling pathways are redundant due to the redundancy of pathway terminal transcription factors (TFs) was tested by screening 35 signalling ligands (SLs) for rescue of a decapentaplegic (dpp) hypomorphic wing growth phenotype. The screen identified three examples of partial rescue: Hedgehog (HH), Semphorin 1a (SEMA1A) and Wnt ortholog 2 (WNT2). HH overexpression with dppGAL4 may increase the expression of DPP activity from the hypomorphic dpp alleles. However, SEMA1A and WNT2 did not phenocopy ectopic expression of HH or DPP and neither SEMA1A nor WNT2 were required for wing growth suggesting substitution of DPP for partial restoration of wing growth. The WNT2 rescue was dependent on the Frizzled 4 (FZ4) WNT receptor excluding the possibility that WNT2 weakly binds the DPP receptor. Although examples of phenotypic nonspecificity of SL function were identified, this is an expectation, and not direct proof, of the hypothesis of TF redundancy. Screen Report SummaryAn expectation of a hypothesis proposing that cell-to-cell signalling pathways are redundant due to the redundancy of the pathway terminal transcription factors was tested by screening for replacement of one signalling ligand (DPP; SLa) with another SLb for wing growth. Three non-DPP SLs were identified in the screen of 35SLs: HH, SEMA1A and WNT2. Genetic analysis of Sema1a and Wnt2 suggests functional complementation of dpp for wing growth suggesting that SEMA1A and WNT2 partially replace DPP for wing growth. Therefore, an expectation of the hypothesis is met.

Cellular metabolism and gene transcription are closely linked. The conserved transcriptional regulator SIN3 acts as a scaffold for histone deacetylase (HDAC)-containing complexes and is crucial for development, stress resistance, and overall organismal health. SIN3 regulates metabolic gene expression in Drosophila cultured cells, however, an understanding of the extent of its role in coordinating responses to metabolic stress in whole organisms is incomplete. In this study, we explored how SIN3 controls glycolytic gene expression across developmental stages and under genetic and dietary disruption of glycolysis in Drosophila melanogaster. Focusing on four key glycolytic enzymes: phosphofructokinase (Pfk), enolase (Eno), pyruvate kinase (Pyk), and pyruvate dehydrogenase beta (Pdhb), we found that reducing Sin3A levels increases their expression in both larvae and adults, indicating that SIN3 plays a consistent role in balancing metabolic gene transcription. Genetic interaction experiments indicate that Sin3A interacts with Pyk and Eno, regulating transcription in a gene-specific manner. Disrupting glycolysis via genetic or dietary means alters glycolytic gene expression, and SIN3 modulates this response. These findings indicate that SIN3 functions as a metabolic sensor, regulating transcription in response to cellular metabolic stress. Additionally, we demonstrate that reducing Sin3A levels shortens Drosophila lifespan on both low- and high-sucrose diets, emphasizing the importance of SIN3 in longevity. Overall, these results show that SIN3 is a context-dependent regulator of glycolytic gene expression and lifespan in Drosophila, integrating metabolic signals with chromatin-based transcriptional regulation. SummaryTo survive and thrive, organisms must adapt to distinct metabolic inputs. We investigated the response of the conserved transcriptional regulator SIN3 to metabolic stress and its control of glycolytic gene expression in Drosophila melanogaster. By measuring glycolytic gene expression, testing genetic interactions, and assessing lifespan under genetic and dietary perturbations, we found that Sin3A knockdown elevates glycolytic gene expression in a gene-specific manner and decreases longevity. SIN3 also modulates transcriptional responses to disrupted glycolysis and influences lifespan under sucrose stress. These findings identify SIN3 as a context-dependent transcription regulator that links gene expression with organismal metabolic adaptation.

Determining spatial and temporal gene expression is crucial for understanding animal development and physiology. Promoter reporters are powerful tools used to dissect how cis-regulatory elements and trans-acting factors control gene expression. Many fluorescent proteins used in promoter reporters, however, have long half-lives (>24 hr) which limit the study of dynamic expression. Destabilizing sequences like PEST reduce the half-life of reporter proteins to provide a more representative readout of gene expression. mlt-11 is a putative protease inhibitor known to oscillate in expression at the mRNA and protein level, yet a mlt-11p::mNeonGreen::3xFLAG::PEST::tbb-2 3UTR promoter reporter did not detectably oscillate. We systematically dissected this transgene, finding that placement of a 3xFLAG tag adjacent to a PEST sequence severely hampered oscillatory expression of a mNeonGreen promoter reporter. Surprisingly, reporter designs that effectively oscillate with a GFP or mNeonGreen fluorophore fail to oscillate when mStayGold is used, with these reporters remaining detectable over 24 hours following promoter inactivation. In addition, other tested epitopes (Myc, ALFA, OLLAS) did not hamper PEST-dependent destabilization but led to varying levels of reporter expression. This study details key considerations for designing destabilized fluorescent promoter reporters.

Pennycress (Thlaspi arvense L.) is an intermediate winter oilseed crop that has only recently been domesticated for agronomic use. Improving agronomic traits requires sources of genetic variation, and mutagenesis is frequently used to help overcome the limitations of natural populations. We investigate the impact of Ethyl methanesulfonate (EMS) on genetically effective cells (GECs) to characterize the intra-individual genetic variation of EMS mutagenesis in pennycress. We identified that pennycress contains at least 4 GECs which, when treated with EMS, create unique mutations across different branches within the same individual plant. We then propagated the M2 plants for whole genome sequencing, providing extensive characterization of the EMS mutation profile and developing a gene index as a resource for future reverse genetic screenings. Article SummaryPennycress is an emerging winter oil seed crop in the American Midwest. Domestication efforts have advanced rapidly through a combination of genetic techniques. One of the most successful methods has been the use of a mutant gene index, a large collection of pennycress seed where new genetic variation has been created through Ethyl methanesulfonate (EMS). EMS mutations are not uniform however, and a single treated seed can have wide genetic variation within the resulting plant. We investigate the role of genetically effective cells on EMS variation, and present the full EMS population as a resource for further pennycress domestication efforts.

O_LIWhen species hybridize, resistance to introgression is presumably due to selection against hybridizing alleles. While many studies have characterized direct selection at these sites, alleles may resist introgression through correlational selection. Here we investigate the role of direct and correlational selection in reducing introgression at the color locus in Ipomoea cordatotriloba. C_LIO_LIWe used recombinant inbred lines that varied in limb color, flower size and sugar concentration to estimate the fitness advantage of the flower color allele via direct and correlational selection. To assess the effect of correlational selection on fitness, we ask if floral size or nectar sugar concentration is correlated with fecundity in pink- but not white-limbed lines. C_LIO_LIWe find no evidence for direct selection on flower color across four fitness components - germination, survival, fecundity, and siring success. Instead, both flower size and sugar concentration significantly correlate with fecundity in pink, but not white limbed lines. As a result, correlational selection on the color allele opposes introgression when recurrent migration is low (<3%). C_LIO_LIThese results demonstrate that correlational, rather than direct, selection is sufficient to resist introgression via hybridization and suggest that correlational selection is an underexplored mechanism to generate resistance to introgression across multiple loci. C_LI

Arabidopsis thaliana naturally occurs across a wide geographic range and displays extensive natural variation in several traits including adaptive responses to the abiotic environment (e.g. temperature, drought, salt). Quantitative techniques like Genome Wide Association Studies (GWAS) enable mapping the genetic basis of such environmental responses and benefits from extensive genetic variation, but the size of the chosen diversity panel is often limited by phenotyping capacity. Most studies therefore use subpanels, often based on maximization of genetic diversity. However, this type of selection may overrepresent cosmopolitan alleles and underrepresent rare environment-specific alleles. Here, we demonstrate that the genetic variation in a GWAS subpanel of Arabidopsis thaliana accessions depends almost entirely on the number of accessions in the panel and very little on the composition of the panel. We present the EcoCore panel designed by grouping accessions of the 1001 genomes (1001G; 1135 accessions) collection, based on their native collection environment and selecting an equal number of accessions from each environment. We assessed hypocotyl lengths of plants grown at control and ambient high temperatures (20{degrees}C and 28{degrees}C) for 913 accessions of the 1001G and mapped these traits with the full 1001G panel versus the EcoCore panel. The EcoCore panel revealed novel genetic associations with hypocotyl length which is attributed to enrichment of alleles from rare environments. We present the EcoCore panel as a manageable resource for studying phenotypic plasticity and the genetic basis of plant-environment interactions.

Crossover recombination events during meiosis repair double-strand DNA breaks and ensure accurate chromosome segregation in most organisms. For many species, the genomic distribution of crossovers is nonrandom and sexually dimorphic. While many species evolved kilobase-scale "hotspots" for crossover formation, the Caenorhabditis elegans genome lacks hotspots, and crossovers are enriched across megabase-scale domains. Further, genetic and cytological studies indicate the crossover frequency in C. elegans spermatogenesis is higher relative to oogenesis in many but not all genetic intervals. To determine the genomic features that contribute to the sexually dimorphic recombination landscape in the absence of hot spots, we defined and analyzed the recombination landscape across the whole genome in C. elegans using whole-genome sequencing and high-resolution recombination mapping in single worms bearing recombinant chromosomes from individual sperm and oocytes. We find that the spatial distribution of crossovers is sexually dimorphic on chromosomes I, II, and III, and that the global rate of double-crossover events is 4.7-fold higher in spermatocytes. Additionally, we find that pairing and synapsis may contribute to the sexually dimorphic crossover landscape. In comparison to the spermatocyte crossover landscape, a higher proportion of oocyte crossovers are formed in the domains directly adjacent to the pairing centers of each chromosome. Further, reducing the genetic dosage of the synaptonemal complex central region protein SYP-2, which is a meiotic chromosome structural protein required for homologous chromosome synapsis, reshapes the oocyte crossover landscape to resemble observations in wild-type spermatocytes. Finally, we found that spermatocyte crossovers are partially enriched in H3K36me3-marked euchromatic regions, while many oocyte crossovers are enriched in H3K27me3-marked heterochromatic regions. Taken together, our studies reveal how synaptonemal complex component dosage and local chromatin states influence crossover placement and the sex-specific regulation of meiotic recombination. Author SummaryProduction of viable eggs and sperm depends on accurate chromosome segregation during meiosis. Segregation of parental copies of homologous chromosomes requires the reciprocal exchange and physical linkage of DNA that arises through crossover recombination. Increasing evidence indicates the existence of sexual dimorphisms during meiotic recombination. In this study, we generated and analyzed high-resolution recombination maps specific to spermatogenesis and oogenesis in the nematode C. elegans, which reveals sex-specific crossover distributions and a higher rate of crossing over in sperm cells. Further, we indicate how specific chromosomal features and structures differentially affect the crossover landscape in eggs versus sperm. Our work highlights how, in a system absent of pre-defined "hotspots" for recombination, local chromatin structures, chromosomal pairing domains, and the abundance of synaptonemal complex proteins are potential drivers for establishing the observable sex differences in crossover recombination.

Two different individuals can consume an identical diet but experience different physiological outcomes. While there are many potential mechanistic reasons for this, one increasingly recognized reason is that differences in an animals microbiome can lead to differences in the processing of dietary nutrients. Thus, though a diet might start out the same, how it is experienced by a host is dependent on their microbiome. While exciting, mechanistic studies of diet-microbiome-host effects are often limited by the lack of high throughput laboratory techniques to identify and define interactions between dietary metabolites, microbial metabolism, and host biology. We hypothesized that the model organism Caenorhabditis elegans is an advantageous animal model for rapidly identifying and genetically dissecting interactions between dietary nutrients, microbial metabolism, and host physiology. Here, we used an established model of the effects of dietary glucose on insulin resistant mutant animals (daf-2/IR mutants) to study how differences in bacterial metabolism influence the consequences of dietary sugars on animal physiology. We found that the effect of dietary sugars on daf-2 mutant physiology is dependent on how the microbiome metabolizes dietary sugars. We found dietary sugar suppresses multiple daf-2 mutant phenotypes in the presence of some bacteria but has no effect in the presence of others. To determine how bacteria mediate the effects of dietary sugars on host physiology we screened 5,000 transposon mutations in the canonical C. elegans dietary bacteria, E. coli OP50 for effects on animal insulin signaling. From this, we found that the effects of exogenous sugars on the phenotype of daf-2 mutant animals is dependent on the function of pyruvate dehydrogenase in bacteria and that the loss of bacterial pyruvate dehydrogenase genes (ex. aceE) is sufficient to mimic the effects of dietary sugars on dauer formation, longevity, and gene expression in insulin signaling deficient animals. Collectively, our findings further support the growing body of evidence that the effects of dietary nutrients on animal physiology can be influenced by the gut microbiome. In addition, these studies demonstrate the advantages of the C. elegans model system for studying 3-way diet-microbiome-host interactions that are difficult to dissect in other model systems.